Identifying a Two-State Hamiltonian in the Presence of Decoherence 
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Mapping the system evolution of a two-state system allows the determination of the effective 
system Hamiltonian directly. We show how this can be achieved even if the system is decohering 
appreciably over the observation time. A method to include various decoherence models is given 
and the limits of this technique are explored. This technique is applicable both to the problem 
of calibrating a control Hamiltonian for quantum computing applications and for precision experi- 
ments in two-state quantum systems. The accuracy of the results obtained with this technique are 
ultimately limited by the validity of the decoherence model used. 
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I. INTRODUCTION 

Recent developments in the fabrication and control of 
few state quantum systems have allowed unprecedented 
tests of both our experimental and theoretical under- 
standing of their behaviour. This is typified by the con- 
tinuing interest in quantum information processing and 
efforts to construct a quantum computer, but equally 
applies to other more mainstream technologies such as 
atomic clocks and quantum optics 0, 0, Q . 

Precision control of these systems inevitably requires 
very accurate information about the system Hamilto- 
nian. Traditionally this has been obtained via a variety 
of experimental and theoretical techniques, though as we 
move toward even higher precision experiments, the need 
for more efficient characterisation techniques is of utmost 
importance. 

Recently, we introduced the concept of characteris- 
ing an unknown two-state quantum system by mapping 
its time evolution and using this data to determine the 
underlying Hamiltonian 0, Mapping the evolution 
in this way has the advantage that minimal knowledge 
of the system is required initially. This is in contrast 
to tomographic methods which are in common use to 
map the fidelity of a given quantum operation, which 
assumes basic knowledge of the system, specifically the 
ability to rotate into different input and measurement 
bases HUSH. 

For solid-state systems, characterisation is especially 
important, as the effective Hamiltonian is strongly de- 
pendant on the fabrication process and control systems. 
The result is that the Hamiltonian of otherwise identical 
systems can vary significantly from device to device and 
even between different regions of the same device, making 
calibration or characterisation of the device a necessity. 



Typically tomography is applied to a system whose be- 
haviour is known, but where the decoherence processes 
need to be characterised. Here we take the opposite ap- 
proach in which the dominant decohering mechanisms 
are known, but it is the strengths of the various deco- 
herence channels that each system experiences, as well 
as the system control parameters which are unknown or 
not known to sufficient precision. In this paper we ex- 
tend our previous results [5| by demonstrating how the 
rate of decoherence can be characterised along with the 
system Hamiltonian. 

Initially, we consider the effect of a pure dephasing de- 
coherence channel (section llllf) and how this can be in- 
cluded in the characterisation process. This simple model 
is then extended to a more general decoherence model 
in section IIVI and we derive analytic expressions for the 
Fourier transform of the time evolution of the system. 
Finally, in section we consider the way in which the 
resulting uncertainties scale as a function of the num- 
ber of measurements and discuss the use of our analytic 
results in fitting experimental data. 



II. SINGLE QUBIT CHARACTERISATION 

Consider a two-state quantum system with a control 
Hamiltonian of the following form 



H = -[sm(8)a x + cos(6)a z ], 



(1) 



using the unnormalised Pauli matrices (a), where the 
z-axis is defined by the measurement basis. The angle 
9 parameterises the ratio of the a x and a z components 
and d is the magnitude of the Hamiltonian. When the 
control Hamiltonian is turned on, the evolution of the 
z-projection of the state of the system which is initially 
in the state z(0) = 1, is given by 
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z{t) = cos(di) sin 2 (6>) + cos 2 (6>). 



(2) 
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We have the freedom to choose the alignment of the 
x-axis and therefore do not include the effects of the 
dy component initially by arbitrarily setting the a;-axis 
to coincide with the Hamiltonian. Once an initial axis 
is defined either by characterisation or on experimen- 
tal grounds, then the angle between any further control 
Hamiltonians and the a;-axis can be determined using ad- 
ditional measurements 



To measure z{t) experimentally, typically the system 
is initialised in some known state, allowed to evolve for 
some time At and then measured. The amount of time 
the system is allowed to evolve for is progressively in- 
creased in increments of At, giving the evolution, z(t), as 
a function of time at N t time points over a total obser- 
vation time i b = N t At. The process must be repeated 
N e times to determine an ensemble average as a function 
of evolution time. This is the conventional coherent os- 
cillation experiment and results in a total of Nt — N e N t 
measurements of the system which will, in general, be 
very large. In this paper we will only consider so-called 
'strong' measurements where the measurement projects 
the system onto one of the two available measurement 
basis states. 

The various parameters of the Hamiltonian can be de- 
termined from the normalised Fourier spectrum of the 
time evolution, .F[z(i)]. As the evolution is purely sinu- 
soidal, in the limit of no decoherence the Fourier spec- 
trum is particularly simple. It comprises ^-functions at 
frequencies lu = and u = d, with magnitudes Fq = F(0) 
and F p = F(d) respectively. An example spectrum is 
shown in Fig. ^ with the appropriate features labelled. 
The position of F p gives the value for d, whereas 9 can 
be found from the peak heights using 



cos 2 (6>) = Fq = 1 



2F p . 



(3) 



If there is an inherent measurement error probability r\ 
(or equivalently some probability of initialisation error), 
this results in a reduction of the amplitude of the oscil- 
lations in a well defined manner. This effect can be com- 
puted by calculating the sum of the Fourier spectrum, 
which for the case of pure delta functions is 



Y,^W)] = F + 2F p = l-2r 1 . 



(4) 



The noise floor apparent in the Fourier spectrum stems 
from both the inherent noise in the experimental setup 
and the discretisation or 'projection' noise ^3 due to 
the fact that each measurement only returns a binary 
result. As the number of ensemble measurements (-/V e ) 
is increased, the noise floor reduces accordingly. We can 
therefore use the distribution of the noise spectrum to 
assign uncertainty estimates to the various parameters 
measured from the Fourier spectrum 
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FIG. 1: (Color online) Example double sided frequency spec- 
trum showing the labelling of peaks for the case of no decoher- 
ence. The zero-frequency peak F(0) = Fn and the oscillation 
frequency peak F(d) — F p while the effect of only taking a fi- 
nite number of measurements is to produce a noise floor. The 
distribution of the noise floor gives an uncertainty estimate 
for the peak heights and therefore the system parameters. 



III. MODELLING THE EFFECT OF 
DECOHERENCE 

To model the effect of decoherence on the characterisa- 
tion process, we use the Lindbladian formalism. The time 
evolution of the system is then gov erned by the Liouville- 
von Neumann equation [3, llll Il2| , 



dp 
~dl 



where p is the density matrix of the system and 
C(p, Li) = L lP L\ - -{L\L h p}, 



(5) 



(6) 



where Li is the Lindbladian operator corresponding to a 
particular decoherence channel. In general a decoherence 
model can include several different Lindbladian opera- 
tors, each of which corresponds to a different decoherence 
mechanism. 

While the Lindblad formalism allows the inclusion of 
general forms of decoherence, it still assumes the Born 
(weak coupling) and Markovian (uncorrelated noise) ap- 
proximations. These approximations are made when 
analysing decoherence for a range of systemsd Q 
ITfij . There are, however, some systems where it is gener- 
ally believed that the Markov approximation is not valid. 
A notable example of this are systems based on super- 
conducting aubitspit ItH ITsl ITflj . where the dominant 
source of noise is thought to be from background charge 
fluctuators and/or I// noise. In this situation, it is often 
difficult or impossible to write down the evolution of the 
system in the form of Eq. JSJ) and therefore the analysis 
of the decoherence needs to be tailored to the particular 
system in auestionfM I2HI22I I2H I2H . For simplicity, we 
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FIG. 2: (Color online) The z-projection of the time evolution 
of a two state system after preparation in the 2 = 1 state, 
for several different dephasing rates. The magnitude of the 
Hamiltonian (d — 1) is kept constant while the angle between 
the components is varied, (a) For 9 = n/2, the oscillations 
decay exponentially, (b) With 9 = 7r/4, the system under- 
goes a different evolution, but stills decays to the mixed state 
z(oo) — 0. Note the difference in scales on the vertical axes 
in (a) and (b). 

will not treat systems of this form but merely point out 
that our analysis is only valid for systems whose evolu- 
tion can be modelled using a master equation of Lindblad 
form. Alternatively, our approach can be used to deter- 
mine an effective phenomenological model for the system 
if Markovian evolution is assumed. 

Initially we will only consider pure dephasing, as this is 
often considered to be the dominant form of decoherence 
for solid-state qubits [Til l25l l26j . We will consider a 
more general model of decoherence in a later section. The 
Lindbladian operator for pure dephasing is given by, 

L z = \ff~ z a z . (7) 

For the case of pure dephasing, we expect an approxi- 
mately exponential decay in the oscillations given by the 
decoherence rate, T z . As the original evolution consists 
of only one oscillation frequency, the exponential decay 
results in a broadening of the peaks in the Fourier spec- 
trum into Lorentzian (Cauchy) distributions. It is in- 
structive to look at how this behaviour is modified when 
9 =/= 7r/2, both in the time and frequency domains, for 
reasons which will become apparent later. 

Solving Eq. JSJ numerically for pure dephasing, the os- 
cillations do decay exponentially when 9 = ir/2. When 
9 7r/2, the effect of dephasing is no longer purely ex- 
ponential decay, but shows a slow decrease towards the 
maximally mixed state z(oo) = 0. Fig.|21shows the evolu- 
tion of the z-projection for (a) 9 = tt/2 and (b) 9 = 7r/4 
for several different dephasing rates. To characterise a 
system which is undergoing decoherence, we need to be 
able to account for these effects. 

Ideally we would like to solve Eq. JSJ for an arbi- 
trary Hamiltonian to obtain the time domain behaviour, 




frequency [units of d] 

FIG. 3: (Color online) The Fourier transform of the system 
evolution given in Fig. [2] in the presence of dephasing, plotted 
for W > 0. The magnitude of the Hamiltonian, d, is kept 
constant and results displayed for (a) 9 — tt/2, and (b) 9 = 
7r/4. The peak is reduced when 9 7^ 7r/2, as expected, and the 
zero-frequency component increases. As the decoherence rate 
increases, the peaks broaden and the frequency shifts slightly. 



though in general it is non-trivial to invert the evolu- 
tion to compute the Hamiltonian parameters. Instead, 
we will look at the system behaviour in the Fourier do- 
main. Fig. Ogives the Fourier transform of the evolution 
shown in Fig. [21 to demonstrate the dependence on both 
the decoherence rate and the angle 9. 

Observing the peak positions under the influence of 
dephasing, we note that the peak position does not move 
appreciably. As the peak is approximately stationary, 
the magnitude of the Hamiltonian vector can still be de- 
termined from its position (at least to first order). As in 
the case for no decoherence, as 9 is varied from tt/2, the 
zero-frequency peak grows as expected. The ratio of the 
area under the two peaks can be used to obtain a first 
order estimate of the value of 9 while the width of the 
peaks are strongly dependant on the dephasing rate. 

In the case of negligible decoherence, the Hamilto- 
nian parameters were determined from the Fourier trans- 
form of the system evolution and the resulting (S-function 
peaks. If we now consider the situation where the deco- 
herence rate is non-negligible, but slower than the system 
oscillations (T z < d), a broadening of the peaks in the 
Fourier spectrum is introduced. Solving Eq. JSJ) analyti- 
cally provides a functional form which can be used to fit 
the experimental data. Starting from Eq. J5J), we gain 
insight by using the conventional variable substitution 

x = (poi + pio)/2 

V = (poi - Pio)/2« (8) 
z = poo - pn, 

to rewrite the matrix equation as a set of three coupled 
differential equations Q|, where pij are the components 
of the density matrix p. For the case of pure dephasing, 
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this gives 

dx(t) 
dt 

dy(t) 
dt 

dz(t) 

dt 



-dcos(0)y(t) -2T z x(t) 

d[sm(6) z{t) - cos(0) x(t)} — 2T Z y(t) (9) 

-dsia(6)y(t). 



The solution of Eqs. © in the time domain is tractable, 
but does not provide any useful insight due to its com- 
plexity. Instead we take the Fourier transform of the 
set of equations which results in a set of algebraic equa- 
tions in terms of X(lu) — T[x(t)], Y(lu) — T[y{t)] and 
Z(u>) = T[z{t)l 

iujX(ui) = dcos(9)Y(u;)-2T z X(uj) 

iuoY{uj) = d[sa\{9)Z{uj)-cos{9)X{uo)]-2Y z Y(uo) 

iwZ(w) = -dsm{9)Y(cu) + Cjr. (10) 

These can be solved to obtain an expression for Z(u), 
where Cjr is a constant of integration arising from z(0) ^ 
0. This gives the solution to the Fourier Transform of z(t) 
as 



C 



■ rf 2 (2r\+i^)sin 2 (0) 



(11) 



where is still unknown in general. 

It is instructive to consider the case where 9 — tt/2, 
in which we find that setting Cjr = 1 artificially and ex- 
panding to first order around u> = d, the real component 
of Eq. ^2 is given by 



Re[Z(w)] 



(12) 



which is a Lorentzian centred about the frequency w = d 
with width given by T z . 

Returning to the general case, if we assume the system 
starts in the state x(— oo) = y{— oo) = z(— oo) = 0, the 
input state z(0) — 1 can be modelled as an impulse at 
t = 0, resulting in the term 2(0)5(0) being added to the 
equation for dz(t)/dt. Transforming to the frequency do- 
main and taking into account the fact we are using the 
discrete Fourier transform gives 



G 



£(0) 

At 



(13) 



as the contribution from this impulse function. 

For the case of pure dephasing, the system must ap- 
proach the perfectly mixed state in the long time limit, 
i.e. x(oo) = y(po) = z(oo) = 0. This means that the 
boundary condition requirements of the Fourier trans- 
form are automatically satisfied. The steady state solu- 
tion 2(00) = also means the issues of frequency reso- 
lution and phase matching, raised in Ref. 0, are not as 
relevant. As long as enough data is gathered that the 



steady state limit is reached, the time evolution data can 
be 'zero-padded' to increase the frequency resolution. If 
the steady state limit is not reached, then the phase dif- 
ference will need to be minimised as in the case of no 
decoherence, though a residual error will still be present 
due to the amplitude mismatch. 

Providing the necessary boundary conditions are met, 
we now have a general procedure for dealing with a sys- 
tem undergoing dephasing. We take the Fourier trans- 
form of the data, as before, and measure the peak po- 
sition to determine d approximately. We then use the 
analytic form given in Eq. I|llf> and perform a nonlinear 
fit on d, T z and 9 to obtain the system parameters. To 
obtain a more accurate fit, the fitting process can be re- 
peated iteratively using the same equations. An initial 
estimate is obtained for each parameter and then the pa- 
rameters are refitted in turn until the estimates converge. 

The parameters d, T z and 9 predominately control the 
peak position, height and width respectively and the in- 
terdependence of the parameters are second order effects. 
We can see this effect in the numerical results shown in 
section ITTT1 and it is this independence of the parameters 
which insures good convergence, as the covariances be- 
tween the parameters are small. The effect of imperfect 
measurement or initialisation can still be characterised 
by computing the sum of the Fourier spectrum as before, 
see appendix lAl 



IV. CHARACTERISATION WITH A MORE 
GENERAL DECOHERENCE MODEL 

To treat more general decoherence, we add terms which 
model spontanteous absorption (r+) and emission (r_) 
(e.g. thermal population transfer) given by Lindbladian 
operators of the form 



l± = Vr^±- (14) 

The decoherence operator becomes 

C{p, U) = C(p, L z ) + C(p, L_) + C(p, L+), (15) 

i 

containing the three forms of decoherence which, in gen- 
eral, will have three different characteristic rates. To 
apply the method given in section ITTll to the more gen- 
eral case requires a variable substitution to match the 
boundary conditions. We define x'(t) = x(t) — x(oo), 
y'(t) = y(t) — 2/(00) and z'(t) = z(t) — z(oo) and then 
solve as before. The initial conditions must also be re- 
defined such that x'(0) — —x(oo), y'(0) = —y(oo) and 
^'(0) = 2(0) — z(oo). Solving in the steady state limit, 
we get 



x(oo) 

2/(00) 

z(oo) 



2 d cos(6>) y(oo) 

4r z + r + + r_ 

K z(oo) 

r+-r_ 

r++r_ +dsm{9)K 



(16) 
(17) 
(18) 
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where 



The resulting solution in the Fourier domain is 



K = 



2d sin(0) (4T Z + T+ 



4d 2 cos 2 (<9) + (4T 2 



(19) 



(O + T+) [1 - z(oo)] - T_ [1 + z(oo)] - d sin(0) [L(w) + L*{-w)\ 



2d? A/ sin 2 (fl) 
M 2 +4d 2 cos 2 (S) 



(20) 



where 



L(w) = 



(Cjf — iw) [y(oo) + iar(oo)] — rfz(oo) sin(0) 



M — lid cos(0) 



m = 2 i lu + 4 r 2 + r+ + r_ 



(21) 



(22) 



and L*(oj) denotes the complex conjugate of L(u>). As 
z(oo) is a constant, Z(uj > 0) = Z'{uj > 0) and Z(0) = 
Z'(0)+z(oo). This solution is an algebraic combination of 
the free variables (d, 9 and T's) and can therefore be used 
as a fitting function for the transform of the oscillation 
data, as outlined previously. 

The effect of different decoherence channels can be dif- 
ficult to discriminate in the Fourier domain, due to their 
similar action on the spectrum. Hence a multi-parameter 
fit will exhibit large covariance terms, with a relatively 
flat potential surface of the fitting function that will not 
easily converge. In this case the underlying physics of 
the system should be used to connect the different deco- 
herence rates and their asymptotic values using thermo- 
dynamic, or other physical arguments. In this way we 
can use the physics to provide additional constraints to 
improve the success of the fitting procedure, an example 
of this is now given for the model in question. 

To illustrate how this process works, we consider the 
special case when the measurement (initialisation) axis is 
coincident with the axis in which the dephasing acts. In 
this situation we can measure the effect of the other (non- 
dephasing) decoherence terms separately and therefore 
reduce the number of free parameters. By repeating the 
experiments as detailed earlier, in the limit of either d = 
or 9 = 0, we build up a picture of the non-Hamiltonian 
evolution. As the system does not have a mechanism 
to move away from the z-axis, the influence of pure de- 
phasing is effectively removed and the population decay 
results purely from the absorption and emission terms 
only. This situation corresponds physically to the limit 
of either no driving field (d = in the rotating frame) or 
the large detuning limit (9 = 0) where the system eigen- 
states are coincident with the measurement basis (<r z ). 
The system evolution in this limit is illustrated in Fig. 0] 
for an example system where T_/r + = 5. Notice that 
the path taken by the system is different depending on 
which state is used for initialisation, though the steady 
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FIG. 4: (Color online) The evolution of a two-state system 
under the influence of both spontaneous absorption and emis- 
sion, in the limit of large detuning. For this example, the 
emission rate is five times that of spontaneous absorption 
(r_/r+ = 5). The path taken by the system depends on 
which initial state is used, though the asymptotic behaviour 
is the same. The two paths are labelled z-i and z\ depending 
on whether the system is initialised in the ground (z(0) = —1) 
or excited state (z(0) = 1) respectively. 



state limit, z(oo), is the same for both. The steady state 
population is given by 



z(po) 



• r. 



(23) 



which provides one equation for determining the two 
rates. Note Eq. (|53J) is just Eq. (JTSJ with d = or = 0. 
This means that by observing the long time behaviour 
only we can reduce the number of free parameters by 
one. We can use this result to define the ratio of T_ to 
r + even in the presence of dephasing and Hamiltonian 
evolution. 

Observing the time behaviour of the total system, 
Eq. (5), with d — or 9 = provides us with another 
handle, as shown in Fig. 0] If we label the trajectories 
taken by the system from the two initial states as Z—\ and 
Z\ for the ground and excited states respectively, we can 
then fit the curves to determine both decay rates. Alter- 
natively, plotting the difference between the trajectories 
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gives a simple expression, 

z 1 (i)-z_i(i) = 2exp[-t(T + +r_)] 3 (24) 

which can be easily fitted to determine the sum of the 
rates. This can then be used with Eq. I|23[) to deter- 
mine both rates independently without the need to fit a 
double exponential. In practise, the experiment could be 
conducted by binning each measurement result based on 
the previous measurement and therefore the initialisation 
state. This would save time in the initialisation phase as 
single qubit operations could then be minimised. 

Using this type of auxiliary experiment, the number 
of free parameters in the expression for Z(ui) can be re- 
duced, resulting in better convergence during the fitting 
process. This in turn gives higher precision estimates for 
the system parameters for the same number of measure- 
ments. 



V. ESTIMATING THE UNCERTAINTY 

To illustrate the method developed so far, we simulate 
an experiment with an arbitrary example Hamiltonian 
H = 0.93(7^ + 0.38cr 2 (d = 1, 8 = 1) and include a 
pure dephasing term with decoherence rate Y/d = 0.1. 
We numerically solve the system evolution for i Q b = 15 
and then include the effect of finite measurement by sim- 
ulating projective measurement, using N t = 1000 and 
N e = 50. The resulting Fourier spectrum is then fitted 
using the Levenberg-Marquardt nonlinear regression al- 
gorithm 27, 28] to perform a non- linear fit in the Fourier 
domain, see Fig. |5] 

Using conventional nonlinear fitting routines to fit the 
measured data with the functions developed in section lTvl 
has the advantage that these routines also provide un- 
certainty estimates based on the goodness-of-fit. These 
uncertainties can be used to directly estimate the un- 
certainty in the final Hamiltonian parameters using the 
same relations derived for no decoherence [5j. We obtain 
rough estimates for the various parameters and then per- 
form an iterative fitting process where each parameter is 
varied in turn. Generally, convergence is achieved within 
2-3 interations, although this depends on the number of 
measurements and ultimately on the quality of the ex- 
perimental data. The numerical values obtained from 
the data shown in Fig. [S] are given in Table [I] for an ex- 
ample run. The true value (x), the estimate (x) and the 
3cr confidence interval (Sx) are given for d, 9 and Y z . 

As in the case of no decoherence, we are ultimately 
interested in how the parameter uncertainties scale with 
increasing number of measurements. Fig. shows the 
fractional uncertainty (la level) for the a z component 
of the Hamiltonian and the decoherence rate of the ex- 
ample system, as a function of the number of measure- 
ments. Similar to the case for no decoherence, the fre- 
quency uncertainty is predominantly controlled by the 
time resolution and is typically much smaller than the 
other uncertainties. Using a constant number of time 
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FIG. 5: (Color online) Simulated data for evolution due to 
the example Hamiltonian, H = 0.93(73, + 0.38cr z . (a) The 
time series data for N e = 50 (points) which approximates the 
ensemble average of the evolution (solid line) . (b) The Fourier 
transform of the time domain data (points) plotted with the 
fitted function (solid line) using the estimated parameters in 
Table □ 
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TABLE I: Example values from a simulated run of the fitting 
procedure discussed above, using the data shown in Fig. 
The true value (x), its estimate (x), the uncertainty (Sx) and 
the fractional uncertainty (5x/x) are given for the three sys- 
tem parameters d, 9 and F z with Nt = 1000 and N e = 50. 



points (Nt = 1 x 10 6 ) and increasing the number of en- 
semble measurements (N e ) we find that the fractional un- 
certainty scales proportional to 1 / yj Nt ■ The fractional 
uncertainty is approximately a factor of 7 larger than 
for no decoherence, which is equal to 1/^/Nt M This 
represents the penalty for fitting including the decoher- 
ence terms. The uncertainty does not change appreciably 
when the decoherence rate changes, with the curves for 
Y/d = 0.1, Y/d = 0.01 and Y/d = 0.001 giving identical 
behaviour. 

The limiting factor in both the fitting procedure and 
the uncertainty analysis is, 'how appropriate is the deco- 
herence model?' If the model used is not suitable, this 
will be apparent as the fitting procedure will not satisfac- 
torily converge, even after many iterations. This can be 
most easily determined by inspection of the Fourier spec- 
trum shape, when compared to that of the closest fitting 
parameters. If the model is found to be the limiting factor 
in the parameter estimation, then a more sophisticated 
model is required which may require additional experi- 
ments. An example of this type of additional experiment 
is the process discussed earlier for determining sponta- 
neous absorption and emission terms independently from 
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FIG. 6: (Color online) The uncertainty estimate as a func- 
tion of number measurements for an example Hamiltonian 
H = 0.93(7^ + 0.38rr 2 undergoing pure dephasing in the bare 
qubit basis. The fractional uncertainty for the la confidence 
interval is plotted for both the a z component of the Hamil- 
tonian and the dephasing rate F z . The scaling is approxi- 
mately proportional to 1/^/Nt and the absolute fractional 
uncertainty is independent of the decoherence rate. 



the effects of dephasing. In situations where the Born 
and/or Markov approximations break-down, it may not 
be possible to model the decoherences using a closed-form 
expression for the time evolution, such as Eq. JSJ • In this 
case, a more sophisticated analysis would be required to 
determine the exact decoherence processes. 



VI. CONCLUSION 

Mapping the time domain evolution of a two state sys- 
tem has previously been shown to provide a systematic 
method for characterising the system Hamiltonian for 
quantum information processing applications. The ma- 
jor drawback being the long decoherence times required 
to accurately map the system evolution. In this paper, 
we have shown how this technique can be extended to in- 
clude the case where the coherent oscillations are damped 
due to the effects of decoherence. 

The effects of relaxation and dephasing are considered 
and a procedure developed to incorporate general deco- 
herence models. This provides an 'after the fact' analysis 
technique which also allows for arbitrary accuracy char- 
acterisation given enough measurement data. The tech- 
nique can be applied to two-state experiments other than 
those required for quantum computing applications and 
the achievable precision is only limited by the validity of 
the decoherence model used. 

The procedure requires deriving a fitting function 
based on a model for the decoherence processes affecting 
the system. We derive the result for a simple model in- 
cluding dephasing and spontaneous absorption and emis- 



sion, acting in the measurement basis of the physical 
qubit. Other more complex models can be included, as 
long as the model can be written down as a closed form 
master equation. In cases whether the decoherence can- 
not be written in this form, the procedure can be used 
with an effective model and the convergence of the pa- 
rameter estimates provides an indication of the validity 
of this effective model. 

The complexity of the decoherence model which can 
be employed is restricted by the need to fit the result- 
ing measurement data. If too many degrees of freedom 
are introduced into the model, the fitting procedure will 
not converge sufficiently to provide usable parameter es- 
timates. In this situation, we find that using the physics 
of the system and performing other auxiliary experiments 
can provide additional constraints on the parameters and 
therefore reduce the uncertainties. This technique is ulti- 
mately limited by the accuracy of the decoherence model, 
the stability of the Hamiltonian and the number of mea- 
surements taken. 
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APPENDIX A: DETERMINATION OF AN 
INITIALISATION OR MEASUREMENT ERROR 

First we show that the initial state of the system, z(0), 
can be determined from the Fourier transform of the 
data. The definition of the discrete Fourier transform 
is given by 



JV-l 

T[z{k)] = z(kAt)e 1 ' 

k=0 



:kn/N 



(Al) 



where N is the number of time or frequency channels, k 
denotes the channel number of the time series, n denotes 
the frequency channel number and both n and k are in- 
tegers. If we compute the sum of the discrete Fourier 
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spectrum over all frequency channels, this gives 

N-l N-l N-l 

n=0 n=0 k=0 

Interchanging the order of the summations gives 

N-l N-l N-l 



(A2) 



/N 



(A3) 



n=0 



fc=0 



n=0 



as z(k At) is independent of n. If we consider the inner 
summation term, evaluating the real and imaginary parts 
separately gives 



N-l 



N-l 



N-l 



e ^kn/N = J2 cos(Lukn/N) + i J2 sin{ujkn/N). 

(A4) 



n=0 n=0 

When k = 0, the cosine term is 
N-l 



n=0 



cos(ujkn/N) 



n=0 



= N, 



(A5) 



fc=0 



whereas, when k > the summation over an entire period 
cancels out for suitable values of N, giving 



N-l 



cos(uikn/N) 



(A6) 



k>0 



Similar, the imaginary sine term is equal to zero for all 
values of k. Putting this together we find that the inner 
summation is equal to 



N-l 



e iukn/N _ NSk Q 



(A7) 



n=0 



i.e. the Kronecker delta function. This then gives us 

N-l N-l 

F[z(k At)] =Nj2z{k At)S k ,o = Nz(0), (A8) 



n=0 



fc=0 



so the initial state of the system is given by the sum over 
the discrete Fourier transform of the evolution, divided 



by the number of time points. If the system is under 
going decoherence, the result is the same as long as the 
boundary conditions are still satisfied. If z(oo) ^ then 
the variable substitution discussed in section HV1 can be 
used and the sum computed over the Fourier transform 
of z'{t). 

If we consider some probability rj that the system is 
initialised in the incorrect state (modelled by a 'bit-flip' 
error immediately after initialisation), (ijj(0)\ijj(0)) = (1— 
77) (0| 0) + 7/(l|l), then the initial state becomes z(0) = 
1 — 1r\. Therefore the sum over all frequencies, 



N-l 



£ .F[2(A;At)] = I-277, 



(A9) 



71=0 



gives the initialisation error. 

For some system evolution given by U = exp(— iHt), 
where H is given by Eq. Q, the evolution of the z- 
projection is given by 



z(t) = Ti[a z Up{0)U% 



(A10) 



where p(0) — |0)(0| is the starting state, z(0) — 1. If 
we model an initialisation error as a bit flip before the 
evolution, this gives the evolution of the 'system+error' 
as 

z c P; r °(t) = Tr[a z Uat P (0)a x tf] = -z(t). (All) 

Similarly we find that applying the bit flip error after the 
evolution gives 

Z v° st (t) = Tr[a z 4Up(0)U^a x ] = -z(t). (A12) 

The two error locations are therefore equivalent for the 
Hamiltonian given in Eq. Q and the resulting system 
evolution with an error probability 77 is given by 

z v (t) = (1 - rj)z{t) + f]z CII (t) = (1 - 2r))z(t), (A13) 

which is the measurement error model used previ- 
ously ■ This means that the sum over the Fourier trans- 
form gives the probability of a measurement or initiali- 
sation error or the cumulative effect if both are present. 
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